Irreversible Opinion Spreading on Scale-Free Networks 
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We study the dynamical and critical behavior of a model for irreversible opinion spreading on 
Barabasi- Albert (BA) scale-free networks by performing extensive Monte Carlo simulations. The 
opinion spreading within an inhomogeneous society is investigated by means of the magnetic Eden 
model, a nonequilibrium kinetic model for the growth of binary mixtures in contact with a thermal 
bath. The deposition dynamics, which is studied as a function of the degree of the occupied sites, 
shows evidence for the leading role played by hubs in the growth process. Systems of finite size grow 
either ordered or disordered, depending on the temperature. By means of standard finite-size scaling 
procedures, the effective order-disorder phase transitions are found to persist in the thermodynamic 
limit. This critical behavior, however, is absent in related equilibrium spin systems such as the Ising 
model on BA scale-free networks, which in the thermodynamic limit only displays a ferromagnetic 
phase. The dependence of these results on the degree exponent is also discussed for the case of 
uncorrelated scale-free networks. 
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I. INTRODUCTION 

Over the last few years, the study of statistical and 
complex systems has proved extremely valuable in pro- 
viding insight into many emerging interdisciplinary fields 
of science [H, HI, [||. In particular, many efforts have 
focused recently on the mathematical modeling of a 
rich variety of social phenomena, such as social influ- 
ence and self-organization, cooperation, opinion forma- 
tion and spreading, evolution of social structures, etc (see 
e.g. 0, 1, 1, 0, SIS M EL El)- In this context, much 
attention has been devoted to the investigation of social 
models defined on complex networks, since their topol- 
ogy reflects some key aspects of social structures, such as 
the small-world effect and the high connectivity of local 
neighborhoods [H Q d B3] . 

The interest in complex networks hugely intensified in 
recent years due to the striking similarity observed in 
the structure of many real networks as different as the 
Internet and the World Wide Web, ecological and food 
webs, power grids and electronic circuits, genome and 
metabolic reactions, collaboration among scientists and 
among Hollywood actors, and many others. 

Empirical observations revealed that the degree dis- 
tributions of these networks, P(k), which measure the 
probability that any randomly chosen node has k links, 
are fat-tailed and typically close to a power-law, i.e. 
P(k) ~ fc~ 7 . Unlike the Poisson degree distribution of 
classical random graphs, power-law distributions lack any 
characteristic scales, hence leading to their designation as 
scale-free networks |14| . 

Following these intriguing observations, Barabasi and 
Albert uncovered the fundamental principles that could 
be responsible for the ubiquitous development of scale- 
free structures in natural and artificial systems (l7j . The 
Barabasi- Albert (BA) model is based on the assumptions 



that (i) real networks evolve gradually by the addition of 
new nodes, and that (ii) the mechanism of linear pref- 
erential attachment governs the probability distribution 
of connections between the newly added nodes and those 
already present in the system. Indeed, this simple dy- 
namical model for the formation of networks was shown 
to generate structures with scale- free degree distributions 
with exponent ^ba = 3, while detailed explanations for 
the processes of growth and preferential attachment were 
provided later for several specific cases [13, [13, [HI • 

Besides the extensive and detailed characterization of 
structural and topological properties of many different 
kinds of complex networks, further investigations have 
focused on the dynamical and critical behavior of mod- 
els defined on them. As part of these efforts, different 
spin models have been studied on scale-free networks 
I20I l2ll . [HI, HE [24 1 , as well as on other classes of 

25, HE Him HlUl. 0ntheone 



_ lia 1201, |21|, 123, |23 
complex networks [lj 
hand, these investigations have been motivated by the 
inherent theoretical interest of exploring further proper- 
ties of many standard models of statistical mechanics. 
On the other hand, spin models defined on complex net- 
work structures can contribute to the understanding of a 
rich diversity of processes in contexts as different as ma- 
terials science, sociology, and biology. For instance, as 
regards the application of Ising-type spin models to the 
study of social phenomena, the spin states may denote 
different opinions or preferences, where the coupling con- 
stant describes the convincing power between interacting 
individuals, which is in competition with the "free will" 
given by the thermal noise Q. Moreover, a magnetic 
field can be used to add a bias that could be interpreted 
as "prejudice" or "stubbornness" [28| . 



Cooperative phenomena arising in models defined on 
complex networks are profoundly affected by the inter- 
play between the model dynamics and the underlying 
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network topology. For instance, simulations of the Ising 
model denned on BA scale-free networks reveal that the 
effective critical temperature scales logarithmically with 
the system size, thus showing that only the ordered fer- 
romagnetic phase is present in the thermodynamic limit 
Similar results are obtained for the Ising model de- 
fined on uncorrelated scale-free networks with exponent 
7 = 3 23]. However, for 2 < 7 < 3 the effective crit- 
ical temperatures are found to scale as a power of the 
system size, while for 7 > 3 the finite-size (pseudo) criti- 
cal temperatures converge to finite values in the thermo- 
dynamic limit, hence indicating the occurrence of true 
ferromagnetic-paramagnetic phase transitions [23| . 

In a related context, noncquilibrium order-disorder 
phase transitions were recently observed in the so-called 
magnetic Eden model (MEM) [3ll [32j growing on small- 
world networks [12| . The MEM is a kinetic growth model 
in which the deposited particles have an intrinsic spin. 
In regular lattices, the MEM's growth process leads to 
Eden-like self-affine growing interfaces and fractal clus- 
ter structures in the bulk [3l|. Furthermore, the MEM 
displays an interesting variety of noncquilibrium phe- 
nomena, such as thermal order-disorder continuous phase 
transitions and spontaneous magnetization reversals, as 
well as morphological, wetting, and corner wetting tran- 
sitions [Hj]. 

The aim of this work is to study the irreversible opinion 
spreading within an inhomogeneous society by character- 
izing the dynamical and critical behavior of the MEM 
growing on scale-free networks. According to the growth 
rules of the MEM, which are given in the next Section, 
the opinion or decision of an individual would be affected 
by those of their acquaintances, but opinion changes 
(analogous to spin flips in an Ising model) would not oc- 
cur. Moreover, the structure of scale-free networks pro- 
vides an interesting setting in which the social roles differ 
significantly from one member to another one. Indeed, 
hubs in the scale- free network represent highly influential 
individuals, whose opinions can potentially affect the de- 
cisions of many other individuals in the society. In order 
to represent the underlying scale-free topology, we will 
mainly focus on the Barabasi- Albert model. However, 
the dependence of the results on the degree exponent 
will also be discussed by considering the MEM growing 
on uncorrelated scale-free networks. 

Although this work is mainly motivated by social phe- 
nomena, a magnetic language will be adopted through- 
out. As commented above, physical concepts such as 
temperature and magnetization, spin growth and clus- 
tering, ferromagnetic-paramagnetic phase transitions, 
etc, can be meaningfully re-interpreted in sociologi- 
cal/sociophysical contexts. 

This paper is organized as follows: in Sec. II, details 
on the model definition and the simulation method are 
given; Sec. Ill is devoted to the presentation and discus- 
sion of the results, and Sec. IV contains the conclusions. 



II. THE MODEL AND THE SIMULATION 
METHOD 

According to the model put forward by Barabasi and 
Albert [13 , the scale-free nature of many real networks 
arises from considering two basic dynamical mechanisms: 
growth and preferential attachment. Starting with a 
small number of nodes, m, a new node with m edges is 
added at every time step and is connected to m different 
nodes already present in the system. The probability n 
for connecting the new node to an arbitrary node of de- 
gree k is given by n = k/K, where K = ^ k is the sum of 
degrees taken over all available nodes in the system. Nu- 
merical simulations and different analytic studies, such as 
the continuum theory, the master equation, and the rate 
equation approaches, show that the BA model asymptot- 
ically leads to scale- free networks with exponent ■jba = 3 

After generating one network realization, a randomly 
chosen up or down spin is deposited on a random site. 
Starting from this seed, the MEM growth takes place by 
adding, one by one, further spins to the immediate neigh- 
borhood (the perimeter) of the growing cluster, where the 
interaction energies between neighboring spins shape the 
growth probability distributions. By analogy to the Ising 
model, the energy E of a configuration of spins is given 
by 

where Si = ±1 indicates the orientation of the spin for 
each occupied site (labeled by the subindex i), J > 
is the ferromagnetic coupling constant between nearest- 
neighbor (NN) spins, and (ij) indicates that the summa- 
tion is taken over all pairs of occupied NN sites (i.e., those 
which are connected within the particular network real- 
ization). As with other spin systems defined on complex 
networks, the magnetic interaction between any pair of 
spins is only present when a network edge connects their 
nodes. 

Setting the Boltzmann constant equal to unity (ks 
1 . the probability for a new spin to be added to the 
perimeter of the existing cluster is defined as propor- 
tional to the Boltzmann factor exp(— AE/T), where T 
is the absolute temperature and AE is the total energy 
change involved. In this work, energy and temperature 
are measured in units of J throughout. At each step, all 
perimeter sites have to be considered and the probabil- 
ities of adding a new (either up or down) spin to each 
site must be evaluated. Using the Monte Carlo simula- 
tion method, all growth probabilities are first computed 
and normalized, and then the growing site and the ori- 
entation of the new spin are both determined by means 
of a pseudorandom number. Although the configuration 
energy of a MEM cluster, given by Eq.((T]), resembles the 
Ising Hamiltonian, it should be noticed that the MEM is 
a nonequilibrium model in which new spins are continu- 
ously added, while older spins remain frozen and are not 




FIG. 1: Average degree of occupied nodes for BA networks 
of size N = 1000, m = 3, and different temperatures, as a 
function of the fraction of deposited particles, n OC c/N. No- 
tice that, assuming a constant deposition rate, the fraction 
of occupied nodes is a measure of time. At low tempera- 
tures, the early stages of the growth process are dominated 
by the hubs, since the MEM's dynamics favors the occupation 
of highly connected nodes within ordered neighborhoods. At 
higher temperatures, however, the phenomenon of preferen- 
tial spin deposition is less significant. 
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FIG. 2: Characterization of the deposition process for BA net- 
works of size N = 10 s , m = 3, and different temperatures, as 
a function of the node degrees: (a) fraction of previously occu- 
pied NN sites relative to the growing node degree, unn /k; (b) 
energy change (with opposite sign) involved in the deposition 
process relative to the growing node degree, —AE/k. 



allowed to flip. The growth of MEM clusters proceeds 
by iterating this deposition process until the network be- 
comes completely filled. 

Notice that, since all normalized growth probabilities 
have to be recalculated at each deposition step, the re- 
sulting update algorithm is quite slow. Ensemble aver- 
ages were calculated by considering typically 10 2 — 10 3 
different (randomly generated) networks and 50 — 100 
different (randomly chosen) seeds for each network con- 
figuration. The ensemble-averaging procedures lead to 
very small statistical errors (typically around or below a 
few percent), which can be safely neglected in the context 
of qualitative discussions. Unless otherwise indicated by 
means of error bars, statistical errors of data displayed 
by symbols are comparable to or smaller than the symbol 
sizes. 

The main part of this paper focuses on BA networks, in 
which different values of the BA parameter in the range 
1 < m < 7 are considered. Moreover, results for uncor- 
rected scale-free networks are also presented for degree 
exponents in the range 2 < 7 < 4. Involving a consid- 
erable computational effort, this work presents extensive 
Monte Carlo simulations for BA and uncorrelated scale- 
free networks with sizes up to N = 2 x 10 4 . 



III. RESULTS AND DISCUSSION 

Let us first focus on the dynamical behavior of MEM 
clusters growing on finite BA scale-free networks. In or- 
der to investigate the process of spin deposition, we com- 
pute the degree of the newly occupied node, k occ , every 
time a new particle is added to the system, as well as the 
number of already occupied NN sites in the neighbor- 
hood of the growing node, unn, and the energy change 
involved in the addition of the new spin, AE. 

Figure 1 shows the average degree of occupied nodes 
relative to the network's mean degree, (k occ )/k, for BA 
networks of size N = 10 3 , m — 3, and different temper- 
atures, as a function of the fraction of occupied nodes, 
n OC c/N. According to the MEM's growth rules, at low 
temperatures the system tends to develop highly ordered 
spin domains. Hence, highly connected nodes have larger 
probabilities to be occupied at early times during the 
growth process. The leading role of hubs is clearly ob- 
served in the low-temperature plots of Figure 1, where 
(k OC c)/k 3> 1 at early stages of the growth process, i.e. 
when the number of occupied nodes is of the order of 
a few percent of the total system size. At higher tem- 
peratures, however, the increased thermal noise tends to 
wash out the phenomenon of preferential spin deposition. 
The growth at later times mainly proceeds by the occu- 
pation of less-than-average connected nodes, which leads 
to roughly temperature-independent values of (k occ ) . 




FIG. 3: Mean number of occupied NN sites and mean energy 
change involved in the addition of a new particle to the sys- 
tem, obtained for BA networks of size N = 10 3 and m — 6. 
The temperature-independent relation (njvjv) = k/2 can be 
explained by means of a simple bond representation that mod- 
els the growth process (see text for more details). 



Figure 2(a) shows the behavior of the average number 
of already occupied NN nodes in the local neighborhood 
of a growing node, as a function of its degree, for BA 
networks of size N = 10 3 , m — 3, and different tem- 
peratures. The fraction of occupied NN nodes, upj^/k, 
decreases monotonically with the node degree k, which is 
due to the fact that, as discussed above, highly connected 
nodes are more likely to be occupied at earlier stages of 
the growth process. Moreover, the distributions become 
flatter for increasingly larger temperatures, as expected. 

In the same vein, Figure 2(b) shows —AE/k, i.e. the 
energy change (with opposite sign) involved in the addi- 
tion of a new spin relative to the degree of the growing 
site, as a function of k. In agreement with the previous 
observations, the plots of —AE/k vs k decrease monoton- 
ically and become flatter for higher temperatures. No- 
tice also that the distributions for the lowest temper- 
ature, T = 0.4, are nearly identical in Figures 2(a) and 
2(b), since at low temperatures most spins grow preferen- 
tially parallel-aligned. Instead, the comparison between 
plots corresponding to higher temperatures shows that 
—AE < njvAr, which clearly reflects the onset of thermal 
disorder. 

Furthermore, it is interesting to note that the energy 
change involved in the occupation of highly connected 
nodes is very low and roughly temperature- independent. 
At low temperatures, the hubs become rapidly filled and 
have only a few neighboring spins. Although the number 
of already deposited neighbors increases for higher tem- 
peratures, the net local magnetization, as measured by 
\AE\, remains small due to the increased thermal noise. 

Figure 3 shows the mean values of tinn and AE av- 
eraged over the full growth process, i.e. until the net- 
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FIG. 4: Thermal dependence of the order parameter for BA 
networks of size N = 10 3 and different values of m, as indi- 
cated. 



work becomes completely filled, as a function of temper- 
ature, for BA networks of size N — 10 3 and m = 3. 
At low temperatures, most spins grow parallel-aligned 
and (—AE) ~ (tinn)- However, at higher tempera- 
tures (—AE) decreases towards the (AE) = infinite- 
temperature asymptotic limit, while, instead, (tinn) re- 
mains constant at the value (u^n) — k/2. 

The latter result can be explained by means of a sim- 
ple bond representation for the MEM's growth process, 
as follows. To each pair of occupied neighboring nodes, 
one can assign a directed bond that points from the ear- 
lier occupied site to the later occupied one. Once the 
network is completely filled, an arbitrary node of degree 
k will have a net bond flux <fi (i.e. the difference be- 
tween outwards- and inwards-directed bonds) given by 
<f> = k — 2tinn- Taking network averages, (<f>) = and 
thus (unn) — k/2, irrespective of the temperature. 

So far, we have addressed some dynamical aspects 
of the MEM growing irreversibly on finite-size BA net- 
works. Now, let us investigate the critical behavior of this 
nonequilibrium system. With this aim, we will adopt an 
appropriate order parameter and focus on its thermal de- 
pendence and the behavior of associated distributions of 
probability. These results will be later extrapolated to 
the thermodynamic limit by means of finite-size scaling 
relations. 

The degree of order in a magnetic system can be nat- 
urally characterized by the total magnetization per site- 
i.e., 

M = ^£^. ( 2 ) 

However, according to standard procedures to avoid spu- 
rious effects arising from the finite size of the simulated 
samples (see e.g. [33|, [H, [35[ ) , it is more appropriate to 
consider instead the absolute value of the total magneti- 




FIG. 5: Normalized probability distribution of the magnetiza- 
tion, P(M), for BA networks of size N = 10 3 , m = 3, and dif- 
ferent temperatures. The sharp peaks arising at M « ±0.95 
for T — 1 have been truncated. 



FIG. 6: Effective transition temperatures for BA networks of 
different size and different values of the parameter m (sym- 
bols). Fits to the data using the finite-size scaling relation, 
Eq. ((31), are also shown (solid lines). 



zation per site, \M\, as the order parameter. 

Figure 4 displays the thermal dependence of (\M\) for 
the MEM growing on BA networks of size N = 10 3 
and different values of the BA parameter to. This figure 
shows, on the one hand, that the effect of increasing to 
at a fixed temperature is that of increasing the net mag- 
netization (and, hence, the order) of the system. Indeed, 
since the mean degree is given by k = m(2 — (m+l)/N) w 
2m, the density of connections grows linearly with to and 
tends to order the system. On the other hand, consid- 
ering the effect of increasing the temperature for a fixed 
value of to, we observe that, at low temperatures, the 
magnetization is close to unity and the system grows or- 
dered, while, at higher temperatures, the disorder sets 
gradually on and the magnetization becomes significantly 
reduced. Indeed, this behavior reflects the occurrence 
of thermally-driven, effective order-disorder phase tran- 
sitions in systems of finite size. 

Further evidence for these pseudophase transitions is 
obtained from the thermal behavior of the magnetic sus- 
ceptibility, which can be calculated from the fluctuations 
of the order parameter, as well as from the heat capacity, 
that can be analogously related to energy fluctuations. 
The thermal distributions of these observables, which are 
not shown here for the sake of space, exhibit peaks that 
are correlated to the effective critical transitions observed 
in the magnetization plots. 

Figure 5 shows the normalized probability distribution 
of the magnetization, P(M), obtained for BA networks 
of size N — I0 3 , m = 3, and different temperatures. As 
expected for thermally disordered systems, at high tem- 
peratures the distribution is peaked at M = 0. However, 
as the temperature is lowered, one observes the onset of 
two maxima located at M = ±M sp (0 < M sp < 1), which 
become sharper and approach M = ±1 as T is gradually 



decreased. The smooth shift of the distribution max- 
ima across T ~ T c , from the high-temperature M = 
to the low-temperature nonzero spontaneous magnetiza- 
tion M = ±M sp , is the hallmark of true thermally-driven 
continuous phase transitions [35j . 

These results provide further evidence on the ordering 
effects induced on dynamical systems by underlying com- 
plex network structures, since the MEM growing in ID 
and 2D regular lattices exhibits only pseudophase tran- 
sitions with effective critical temperatures Tf;*f(N) that 
vanish in the N oo limit [3lj . Moreover, the interplay 
between the topology of BA scale-free networks and the 
MEM nonequilibrium dynamics leads to critical phenom- 
ena which are absent in related equilibrium spin systems. 
For instance, as commented above, the Ising model de- 
fined on BA scale-free structures shows effective critical 
temperatures that grow logarithmically with the system 
size, hence displaying only a ferromagnetic phase in the 
thermodynamic limit @, [2(| ■ 

The existence of true critical phase transitions, charac- 
terized by m— dependent finite critical temperatures that 
separate the low-temperature ferromagnetic phase from 
the high-temperature paramagnetic phase, can be fur- 
ther confirmed by means of a standard finite-size scaling 
procedure. 

According to the finite-size scaling theory, developed 
for the treatment of finite-size effects at criticality and 
under equilibrium conditions [36l . \3l\ , the difference be- 
tween the true critical temperature, T c , and an effective 
pseudocritical one, T^^(N), is given by 

\T c -T^f(N)\ ocA^ 1 /", (3) 

where v is the exponent that characterizes the divergence 
of the correlation length at criticality. The effective pseu- 
docritical temperature, T^' (N), is here defined as the 
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FIG. 7: Critical temperature as a function of m (symbols), 
which can be approximated by a linear fit (solid line). The 
inset shows the dependence with m of the exponent v, where 
the dashed line is a guide to the eye. 



value corresponding to (\M\) — 0.5 for a finite system of 
N nodes. 

The symbols in Figure 6 show the effective pseudo- 
critical temperatures for the MEM growing on BA net- 
works of different size and different values of the param- 
eter to. The solid lines are least-squares fits to the data 
obtained by means of the finite-size scaling relation, Eq. 
([3]). The effective pseudocritical temperatures decrease 
monotonically with the network size and lead to finite 
extrapolations in the thermodynamic limit, thus confirm- 
ing the anticipated critical behavior. Notice also that, for 
a fixed system size, the pseudocritical temperatures in- 
crease monotonically with to, as expected. 

Figure 7 shows the critical temperatures obtained for 
different values of to. The monotonic increase of T c (m) 
can be well approximated by the linear relation T c (m) = 
0.522(3) x to + 0.21(1), which is displayed in the Figure 
by a solid line. Moreover, the dependence with to of 
the exponent v is shown in the inset, revealing that v 
decreases monotonically and tends to v s=s 2 for to > 3. 

Let us now consider the irreversible growth of MEM 
clusters on uncorrelated scale-free networks. Since this 
type of networks allows a free choice of the degree expo- 
nent 7, we will investigate how the dynamical and critical 
behavior of the MEM depends on the steepness of the de- 
gree distribution. Moreover, we will compare the critical 
phenomena observed here to previous investigations in a 
related context, namely the Ising model on uncorrelated 
scale-free networks [U, H3] . 

In order to generate uncorrelated scale-free networks, 
we assign to each node a given number k of "stubs" (link 
ends) and match them pairwise at random, with the con- 
straint of avoiding self- and multiple connections. Besides 
the exponent 7, the degree distribution of a finite scale- 
free network is defined by the minimum degree fco and 
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FIG. 8: Average degree of occupied nodes for uncorrelated 
scale-free networks of size N « 10 3 , minimum degree ko = 
3, and different values of temperature and degree exponent, 
as a function of the fraction of deposited particles, n occ /N . 
The phenomenon of preferential spin deposition, which reveals 
that hubs dominate the early stages of the growth process, is 
only significant at low temperatures (as observed in the case 
of BA networks, see Figure 1) and for low degree exponents 
7 < 3. 



the upper cutoff k cu t- Hence, the number of sites with 
degree k is given by N k = (k/k cut )~~< for k < k < k cut , 
and Nk = otherwise, while the total number of nodes 
is N — ^2 k Nk 23]. This procedure, which is based on 
the so-called configuration model (see e.g. [38|, [39|, 140 ]) 
but includes an additional restriction on the maximum 
possible degree of the vertices, was shown to generate 
scale-free networks with no two- and three-vertex corre- 
lations nnm. 

As discussed above, the process of spin deposition dur- 
ing the growth of MEM clusters can be studied by com- 
puting the degree of the newly occupied node, k occ , every 
time a new particle is added to the system (recall Fig- 
ure 1). Figure 8 shows plots of (k occ )/k vs n occ /N for 
uncorrelated scale-free networks of size N w 10 3 , mini- 
mum degree kg = 3, and different values of 7, as indi- 
cated. For each value of the degree exponent considered, 
we calculated the average degree of occupied nodes for 
two different temperatures that characterize the typical 
low- and high-temperature behavior, namely T = 0.4 and 
T = 4. Notice that, due to the procedure followed here to 
generate uncorrelated scale-free networks, the size of the 
networks is controlled by the degree cutoff k cut . Indeed, 
the actual size of the networks represented in Figure 8 
fluctuates by ~ 5% around N = 10 3 . 

The plots in Figure 8 confirm that the phenomenon 
of preferential spin deposition, due to the dominant role 
played by hubs during the early stages of the growth pro- 
cess, is only significant at low temperatures. Moreover, 
this phenomenon is observed to be relevant only for low 
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FIG. 9: Thermal dependence of the order parameter for un- 
corrected scale- free networks of size N w 10 3 , minimum de- 
gree ko = 3, and different values of degree exponent, as indi- 
cated. 



FIG. 10: Effective transition temperatures for uncorrelated 
scale-free networks of minimum degree ko = 3, different size, 
and different values of degree exponent (symbols) . Fits to the 
data using the scaling relations given by Eqs. © and @ are 
also shown (solid lines). 



degree exponents 7 < 3. Since the preferential spin de- 
position is a feature associated to the formation of large 
ordered clusters during the growth process, these results 
agree well, at a qualitative level, with the analogous be- 
havior reported for the Ising model on uncorrelated scale- 
free networks (see e.g. |23j), m which the disorder was 
observed to grow monotonically with the exponent 7. 

Figure 9 displays the thermal dependence of the order 
parameter, (|M|), for uncorrelated scale-free networks of 
size N « 10 3 , minimum degree kg = 3, and different val- 
ues of the degree exponent. As discussed above, the de- 
gree of disorder in the system is observed to grow mono- 
tonically with 7, irrespective of the temperature. More- 
over, for any fixed value of 7, the temperature is shown 
to drive the system through an effective order-disorder 
pseudophase transition. 

In order to further investigate the scaling properties of 
the observed pseudocritical behavior, we can determine 
the effective pseudocritical temperatures for systems of 
different size and extrapolate the results to the thermo- 
dynamic limit. Analogously to the procedure followed 
above (see Figures 6 and 7), symbols in Figure 10 repre- 
sent the effective transition temperatures for the MEM 
growing on uncorrelated scale-free networks of minimum 
degree fco = 3, different size, and different values of the 
degree exponent. The critical behavior turns out to be 
crucially dependent on the steepness of the degree distri- 
bution. Indeed, while the plots with 7 > 3 look similar 
to those obtained previously for the MEM growing on 
BA networks, the plots corresponding to smaller values 
of the degree exponent are observed to diverge, hence im- 
plying the absence of paramagnetic-ferromagnetic phase 
transitions in the thermodynamic limit. 

Using the finite-size scaling relation given by Eq. l[5]). 
the plots for 7 > 3 can be well fitted (solid lines in Figure 



10). For 7 = 3, we obtain T c = 1.14(2) and v = 0.9(1), 
while, for 7 = 4, the scaling relation yields T c — 0.74(1) 
and v = 2.6(2). 

In the same vein, the divergent plots for 7 < 3 can 
be well fitted by a scaling relation analogous to Eq. ([3]), 
namely 

T c e// (A) = A + B x N z , (4) 

as shown by the corresponding solid lines in Figure 10. 
The 7 = 2 plot diverges with exponent z = 0.29(3), while, 
for 7 = 2.5, one obtains z — 0.33(3). 

Comparing these observations with the results re- 
ported for the Ising model on uncorrelated scale-free net- 
works [2l|, HH, we find a qualitative agreement for 
both 7 < 3 (i.e. the absence of a paramagnetic phase in 
the thermodynamic limit) and 7 > 3 (i.e. the existence of 
a finite critical temperature delimiting the paramagnetic- 
ferromagnetic phase transition). Instead, the conver- 
gence of the critical temperature for the 7 = 3 case is 
clearly in contrast with the logarithmic divergence ob- 
served in the Ising model. Indeed, this was also pointed 
out above when comparing the behavior of both models 
defined on BA structures. We can thus stress the fact 
that, for both the MEM and the Ising model, the critical 
features remain unchanged when passing from the BA 
network topology to the uncorrelated scale-free network 
structure with degree exponent 7 = 3. 



IV. CONCLUSIONS 

Dynamical and critical properties of the magnetic Eden 
model growing on Barabasi-Albert scale-free networks 
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were studied by means of extensive Monte Carlo sim- 
ulations. According to former studies on spin models 
defined on complex networks, magnetic concepts such as 
e.g. temperature, spin states, and ferromagnetic cou- 
plings, can be meaningfully reinterpreted in social con- 
texts. Hence, the MEM growing on scale-free networks 
is here proposed as a model for the irreversible opinion 
spreading within an inhomogeneous society, in which dif- 
ferent individuals play different social roles. 

The deposition/spreading dynamics was studied as a 
function of the degree of the occupied sites and showed 
evidence for the leading role played by hubs in the growth 
process. Indeed, it was observed a phenomenon of prefer- 
ential spin deposition, in which highly connected nodes 
tend to dominate the growth dynamics at early times. 
This phenomenon is more significant at low tempera- 
tures, since the MEM's growth rules favor the occupation 
of highly connected nodes within ordered local neighbor- 
hoods. 

Adopting the mean absolute magnetization per spin 
as the order parameter, we observed the occurrence of 
effective pseudocritical order-disorder phase transitions. 
The thermal dependence of the order parameter dis- 
tribution functions showed the characteristic signatures 
of true thermally-driven continuous phase transitions, 
which were indeed confirmed by means of standard finite- 
size scaling procedures. 

The critical behavior found in this noncquilibrium spin 
model contrasts remarkably with the behavior observed 



in related equilibrium spin systems such as the Ising 
model on BA scale-free networks, which in the thermo- 
dynamic limit only displays a ferromagnetic phase. 

Finally, the irreversible growth of MEM clusters was 
studied on uncorrelated scale-free networks. Many of the 
qualitative features in the growth of finite samples were 
found to be common to both BA and uncorrelated scale- 
free networks. However, the latter were observed to ex- 
hibit a very strong dependence on the degree exponent 
7, which manifests itself most significantly in the criti- 
cal properties of the system. The critical behavior for 
7 > 3 was found analogous to previous results for the 
MEM growing on BA networks, while, for smaller values 
of the degree exponent, only the ferromagnetic phase was 
observed to persist in the thermodynamic limit. 

The present findings will thus hopefully stimulate and 
contribute to the growing interdisciplinary efforts in the 
fields of sociophysics, complex networks, and noncquilib- 
rium statistical physics. 
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